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I. INTRODUCTION 



In the recent years, the AdS/CFT correspondence pj-yj has been applied to study strongly 
coupled phenomena in condensed matter physics. Inspired by the idea of spontaneous sym- 
metry breaking in the presence of horizon [4, 5], the Holographic superconductors established 
in 0, 13] are remarkable examples where the Gauge/Gravity duality plays an important 
role. Holographic superconductors were also studied in various backgrounds such as Gauss- 
Bonnet |s- 11], Born-Infeld 12], Horava-Lifshitz theory flfl- 

In the present paper we would like to study a holographic model of superconductor with 

3- 



external magnetic field. 
According to the result of 



15 



'his problem previously has been studied numerically 
15], there exists a critical value of magnetic field, below which a 



charged condensate can form via a second order phase transition. Recently some analytical 
approaches have been proposed to address the universal properties of second order phase 



transitions in holographic superconductors 



17H23I] . In particular, the authors of [18| used 



the variational method for the Sturm-Liouville eigenvalue problem to analytically calculate 
some properties of the holographic superconductors in a (2 + l)-dimensional boundary field 
theory. Here we want to study the effect of external magnetic filed on the holographic 
superconductor analytically by the variational method for the Sturm-Liouville eigenvalue 
problem. To implement a magnetic filed at finite temperature, we consider a magnetically 
charged black hole in AdS^. Then, to incorporate an effectic theory of superconductor, we 
will probe the background with an electric field and a scalar hair. Earlier study showed 
that a holographic superconductor persists in its superconducting state at the temperature 
well below its critical value T c and magnetic field below its critical H c [15j . In this paper, 
we are able to reproduce those critical points using the analytical method mentioned above. 
The paper is organized as follows: in the section II, we recall the setup of Holographic 
superconductor in the external magnetic field. We then develope the analytical method for 
nonzero magnetic field in the section III and study the ciritcal exponent near phase transition 
in the section IV. We summarize our results for various dimensions of condensation at finite 
temperature in phase diagrams in the section V. We also compute critical charge density at 
zero temperature in the section VI. At last we conclude our results in the section VII. 
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II. A HOLOGRAPHIC MODEL OF SUPERCONDUCTOR WITH EXTERNAL 
MAGNETIC FIELD 



In this section, we review the construction of a holographic superconductor in the external 



background 



151 ] . The metric for charged black hole is given by 



dv^ 1 

ds 2 = -f(r)dt 2 + + r 2 (dx 2 + dy 2 ), (1) 
fir) 

where 

*, s r2 M H 2 

f^) = T 2-~ + ^- (2) 
Here black mass and magnetic charge are M and H . The metric is asymptotically AdS± 
with curvature radius L, which will be set to unity without loss of generality. We then 
interpretate H as the external magnetic field in a holographic superconductor living on the 
boundary of the geoemtry. While one of the roots to the equation f(r) = is identified as 
the horizon, r = r + , we can replace M and rewrite 

r 3 tt2 tt2 

f(r)=r 2 -^±- — + \ (3) 

Since the boundary field theory reaches thermal equilibrium with the gravity bulk for a 
static solution, we can identify the Hawking temeprature of charged black hole as the same 
temperature in the dual field theory, that is 

T = . (4) 

An 

To incorporate an effective theory of superconductor, we will porbe the background with 
an electric field and a scalar hair. The asymptotic behaviour of electric field gives chemical 
potential and charge density of condensation, while the vev of scalar hair will tell us whether 
we are in the normal or superconducting phase. The equations of motion were given 



f 2 $ 2 m 2 
2 m 2 ^ 2 

$" + £ $' + HL?_$ = o. (5) 

r f 

We remark that these equations have been solved numerically in the presence of magnetic 
field by the shooting method for the choice m 2 = —2^^. 



3 



III. ANALYTICAL COMPUTATION OF CRITICAL MAGNETIC FIELD 



In this section, we would like to formulate the analytical method in order to study critical 
points. It is convenient to introduce a new cooridnate 

z= r -± (6) 
r 

such that the horizon locates at z = 1 and boundary at z = 0. In this new coordinate, the 
equations of motion of ^ and $ become 

*"M + m = 0, (8) 

where 

f( z ) =r 2 + (z' 2 -h 2 z- z + h 2 z 2 ). (9) 
Here we have introduced a useful normalization of magnetic field, 

v = f£. no) 

To obtain analytical solutions at critical points, we use the ansatze 

P 



$(z)=\r+ n a-z), A 

" +c 



2 > 



*w = Ip^w ( n ) 



where A is one of roots 



A ± = 3 -±J 9 - + m*. (12) 



To bring to the Sturm-Liouville form, we substitute f lTTj) into (J7J) and multiply it with 
function P/j, such that 



with 



[P h {z)F' h {z)]' - + A^ ft (s)F h (z) = 0, (13) 



P h {z) = z 2A f(z), 

Q h {z) = — A(A - l)^ 2A - 2 /(^) " A* 2 *" 1 /'^) - m 2 r 2 z 2A -\ 



/(*) 



Using the Sturm-Liouville method, we manage to minimize the eigenvalue A^ by variation 
of the following expression 

A2 = Jo P h {z)F' h {zfdz + j; Q h (z)F h (z) 2 dz 
h j l Q W h {z)F h {zfdz 

with a trial function Fh(z) = 1 — az 2 subject to the boundary condition F^(0) = and 
normalization Fh(0) = 1. The expression (Ti~4"l) might not have a closed form but can be 
expanded in h 2 order by order as follows: 

, 2 = N (a) + N 2 (a)h 2 

h D (a) + D 2 (a)h 2 + --- + D 2k (a)h 2k + --- 1 ' 

with 

iVo(a) = -(3 — 3a + 5a 2 ), (16) 
6 

N 2 (a) = —(-35 + 63a - 45a 2 ), 
v ; 120 V ; ' 

1 7T 13 

D (a) = -(v^vr - 3 log 3) + (-3 + -= + log3)a + (-— + log3)a 2 , 
6 V3 12 

1 13 343 4 

A(a) = ^(-24 -V3n + 271og3) + (— - v^vr + log3)a + (— - ^= - 31og3)a 2 , 



where D 2k (a) are quadratic function of a for integer fc. If we set /i = 0, we recover the 



analytical result in the s-wave holographic superconndcutor 19fl . One can also obtain more 
precise results by demanding expansion of two or more terms in the trial function, say 
Fh(z) = 1 — az 2 — j3z 3 + 0(z 4 ), then we have to extend 

Ar , , Ar , . 1 a 5a 2 2/3 lla/3 4/3 2 
N (a)^N (a^) = --- + — -- + — + — , 

iV 2 (a)^iV 2 (a,/ 3 )^-i + ^-^ + ^-^-^, 
K J v 1,1 6 10 14 15 7 9 



D k (a) ->• D k (a,f3) 

(17) 

where we have skipped writing down complicated expressions for D k . Now we can view 
that turning on nonzero magnetic field h ^ deforms the eigenvalues Xh away from Ao and 
therefore changes the critical temperature. In the future, without confusion, we will simply 
use A to denote eigenvalues in both cases of zero and nonzero h. 
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FIG. 1: Coefficient 7^ against various h for given m? = —2, A = 2. 



IV. CRITICAL EXPONENT NEAR CRITICAL TEMPERATURE 



Since the condensate (O) is very small near the critical temperature, one could expand 
$(z) as 



*£U A(1 _, ) + g xW . 

Plug into equation (JSJ and integrate by part, one obtains 

X (0) = ArrT I dz- 



(18) 



l-/i 2 ^ 3 -z 3 + /i 2 z 4 ' 

Expand equation ([TBI near 2 = and collect 0(2) term, we have the following relation 



(19) 



r 



T = A - ^2A X (0) 



So we can deduce that 



with 



(0)=l h T c A Jl- 



T 



(20) 



(21) 



7ft 
B-- 



8 / 4tt A 
— m 2 |£> 3 — /z 2 

1 z 2A-2 F(z)2(l _ ^ 



(22) 



,„ ' 1 - h 2 z 3 - z 3 + h 2 z 4 ' 
We obtain 7^ ~ 131 in the limit h — 0, to be compared with the numerical result 144 
obtained in [6J. In the figure 1, we plot coefficient 7^ against various h in the case of 



m 



-2, A = 2 and find that it increases with h. 



V. CRITICAL MAGNETIC FIELD AT FINITE TEMPERATURE 



Earlier study showed that a holographic superconductor persists in its superconducting 
state at the temperature well below its critical value T c and magnetic field below its crit- 
ical H c [15|. In this paper, we are able to reproduce the result using the analytic method 
introduced above. The figure 2 shows that a curve denoting pairs of (T C ,H C ) seperates the 
superconductnig state from normal state. The curve can now be analytically determined by 
(fl4|) for various h. To be explicit, the critical values are obtained via 

rp 3 h\ 

4tt +v 3 ; ' 
H c = r 2 +c h, (23) 

once the eigenvalues are found. 

The two parameters (a, 0) minimization mainl y co rrects the tail part of curve (large h) 
and agrees better with the earlier numerical resultsjlo)] . The decreasing of T c with increasing 
H implies a similar effect to Meissner effect in usual superconductors. This is also a general 
feature for various dimensions of condensate given choices of scalar mass. 



VI. ZERO TEMPERATURE LIMIT AND CRITICAL CHEMICAL POTENTIAL 

A special case occurs for h 2 = 3, where function f(z) has double roots and the charged 
black hole becomes extremal. The horizon of finite size indicates a degenerate ground state 
with finite entropy. Unlike its finite temperature counterpart, we can now have a phase 
transition triggered by quantum fluctuation as charge density varies across some critical 
value. The critical density p c is obtained from ffl4"|) while h 2 = 3. For example, the choice 
of m 2 = —2 allows condensate of two possible dimensions: 



A. A = 1 



In this case, we have the expression 

0.4a + 0.190476a 2 
~ 0.43521 - 0.171785a + 0.0436673a 2 ^ ' 

and the minimum A 2 = —0.33032 with respect to a = —0.837617. This gives the critical 
p c = 0.574735 



m 2 =-2, A=l 



m 2 =-2, A=2 
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FIG. 2: The phase diagram of Tc against He for various A and m 2 . The superconducting phase 
exists in the shaded part below the curve, while normal phase in the upper right part above the 
curve. The dashed (solid) curves is based on minimization of single (two) parameter(s) at constant 
p. We remark that T c has been normalized to 1 when there is no magnetic field in presence. 

B. A = 2 



In this case, we have the expression 

0.4 - 0.190476a + 0.133333a 2 



A 



(25) 



0.0858923 - 0.0873346a + 0.0289007a 2 
and the minimum A 2 = 3.48881 with respect to a = —1.75696. This gives the critical 



Pc 



1.86784 



We also plot critical charge density p c against various upper dimensions A + in the figure 
3 and a linear relation p c oc A + is found for integer and half integer A. The coefficients of 
proportionality (the slope for consecutive data points in the figure 2) are 1.18938, 1.21353, 
1.23341 for one, two, and three parameters minimization respectively. We remark that a 
similar analysis of critical magnetic field for various scalar mass and charge coupling has 



been discussed in 



24]. 



Pc 
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FIG. 3: Critical charge density p c against various dimensions A + for given m 2 . The three pa- 
rameters minimization (blue diamond) gives slightly smaller critical values in comparison to those 
obtained in the two parameters (red square) and one parameter minimization (black circle). We 
remark that the three parameters minimization is obtained by minimizing (a, (3, 7) given the trial 
function F(z) = 1 — az 2 — (3z 3 — 7Z 4 . 

VII. DISCUSSION 

In this paper, we have used the Sturm-Liouville method to analytically compute the 
critical magnetic field for the Holographic superconductor in the superconducting phase. 
At first, the critical exponent is shown to remain 1/2 regardless of external magnetic field, 
as expected from an effective field theory of second order phase transition. The coefficient 
of proportionality 7^ is found to increase with external magnetic field. The phase diagram 
shown in the figure 2 implies an effect similar to the Meissner effect in usual superconductors. 



In addition to reproduce the same phase diagram previously numerically obtained in [15 1 
for the case m 2 = —2, we have found that the same feature persists in cases of other 
permissible m? and therefore condensate of various dimensions. The degenerated ground 
state of Holographic superconductors can be reached by tuning the magnetic field to a specific 
value, where the corresponding charged black hole achieves extremality. One expects a phase 
transition may occur at some critical charge density thanks to the quantum fluctuation. 
Using the analytical method we have found this critical point for condensate of various 
dimensions. We also find that critical value increases linearly with the dimension and the 
coefficient of proportionality is also computed. The Sturm-Liouville method has proved 



useful in analytically study of Holographic superconductors at critical points. We expect its 
application to other critical phenomena which can be stuidied by the holographic method. 

Acknowledgments 

WYW would like to thank the hospitality of Caltech High Energy Theory group in the 
early stage of this project. This work is supported in part by the Taiwan's National Science 
Concil and National Center for Theoretical Science. 



[1] J. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). 

[2] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998); 

|hep-th/9802109[ 
[3] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998). 
[4] S. S. Gubser, Class. Quant. Grav. 22, 5121 (2005). 
[5] S. S. Gubser, Phys. Rev. D 78, 065034 (2008). 

[6] S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, Phys. Rev. Lett. 101, 031601 (2008). 
[7] G. T. Horowitz and M. M. Roberts, Phys. Rev. D 78, 126008 (2008). 
[8] R. Gregory, S. Kanno, and J. Soda, J. High Energy Phys. 0910, 010 (2009). 
[9] L. Barclay, R. Gregory, S. Kanno, and P. Sutcliffe, J. High Energy Phys. 1012, 029 (2010); 
larXiv:1009.199^ hep-th] . 
[10] Q. Y. Pan and B. Wang, Phys. Lett. B 693, 159 (2010). 

[11] Rong-Gen Cai, Zhang- Yu Nie, and Hai-Qing Zhang, Phys. Rev. D 82, 066007 (2010); 

larXiv: 1007.3321 [hep-th]. 
[12] Jiliang Jing and Songbai Chen, Phys. Lett. B 686, 68 (2010) 

[13] R. G. Cai and H. Q. Zhang, Phys. Rev. D 81, 066003 (2010) |arXiv:091 1.48671 [hep-th]]. 
[14] D. Momeni, M. R. Setare, N. Majd, JHEP 1105, 118 (2011). |arXiv: 1003 .03761 [hep-th]]. 
[15] E. Nakano and Wen-Yu Wen, Phys. Rev. D 78, 046004 (2008). 
[16] T. Albash and C. V. Johnson, JHEP 0809, 121 (2008) [arXiv:0804.3466l [hep-th]]. 
[17] X. H. Ge, B. Wang, S. F. Wu and G. H. Yang, JHEP 1008, 108 (2010) [ arXiv: 1002.4901 1 
[hep-th]]. 



10 



[18] G. Siopsis and J. Therrien, JHEP 1005, 013 (2010) [arXiv: 1003.42751 [hep-th]]. 
[19] H. B. Zeng, X. Gao, Y. Jiang and H. S. Zong, JHEP 1105, 002 (2011) [ arXiv: 1012.55641 
[hep-th]]. 

[20] H. F. Li, R. G. Cai and H. Q. Zhang, JHEP 1104, 028 (2011) [arXiv:1103.2833l [hep-th]]. 
[21] R. G. Cai, H. F. Li and H. Q. Zhang, Phys. Rev. D 83, 126007 (2011) [ arXiv:1103.5568l 
[hep-th]]. 

[22] C. M. Chen and M. F. Wu, arXiv:1103.5l30l [hep-th]. 
[23] X. H. Ge. larXiv: 1105.4333 1 [hep-th], 

[24] N. Iqbal, H. Liu, M. Mezei and Q. Si, Phys. Rev. D 82, 045002 (2010) [ arXiv:1003.00T0l 
[hep-th]]. 



11 



